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We classify, up to local unitary equivalence, local unitary stabilizer Lie algebras for symmetric 
mixed states of n qubits into six classes. These include the stabilizer types of the Werner states, 
the GHZ state and its generalizations, and Dicke states. For all but the zero algebra, we classify 
entanglement types (local unitary equivalence classes) of symmetric mixed states that have those 
stabilizers. We make use of the identification of symmetric density matrices with polynomials in 
three variables with real coefficients and apply the representation theory of 50(3) on this space of 
polynomials. 
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I. INTRODUCTION 



Quantum information seeks to exploit quantum states 
as resources for computational tasks, where the phe- 
nomenon of entanglement seems to play a central role in 
quantum advantages over classical protocols. A problem 
at the heart of the study of entanglement is to under- 
stand the essential differences between local and global 
unitary operations on composite systems. Questions such 
as the following arise naturally. How can one characterize 
the nonlocal properties of a given state? Given an input 
state, what is a reasonable description of the set of states 
to which it can be converted using only local operations? 
These questions are known to be hard. Nonlocal prop- 
erties are poorly understood [l|. Regarding the second 
question, the number of polynomial invariants needed to 
determine local unitary equivalence classes of multiqubit 
states grows exponentially with the number of qubits [2| . 

The considerations of the previous paragraph motivate 
the search for entanglement measures and local unitary 
invariants that, being more tractable than a complete set 
of polynomial invariants, will provide a coarser grained 
classification, yet still fine enough to distinguish useful 
states and good entanglement properties. We pursue this 
philosophy with the following scheme. Given an n-qubit 
density matrix p, its local unitary stabilizer subgroup is 
defined to be all the local transformations of state space 
that leave p invariant. If p, p' are local unitary equiva- 
lent, say by a local unitary operation U, then their sta- 
bilizers are isomorphic via conjugation by C/ [3|. Thus 
the conjugacy class of the stabilizer subgroup of a state 
is a local unitary invariant of that state. This invari- 
ant is practical and has proven to distinguish states with 
known useful entanglement properties. Tractability is 
achieved by virtue of the Lie group structure of the local 
unitary group and stabilizer subgroups. Lie groups admit 
analysis through their Lie algebras, which are their tan- 
gent spaces of infinitesimal local operations. In previous 



work, we have classified conjugacy classes of stabilizers 
and local unitary classes of states including the singlet, 
GHZ and its generalizations, symmetric Dicke states (in- 
cluding_the W state and its generalizations) , and Werner 
states 



In this article, we consider the class of symmetric 
mixed states of n-qubit systems. Symmetric states, that 
is, states that are invariant under permutation of subsys- 
tems, are the subject of recent work including: geomet- 
ric measure of entanglement J9l4l3j , efficient tomography 
[14| . classification of states equivalent under stochastic 
local operations and classical communication (SLOCC) 
[l5l4l7| , and our own work on classification of pure sym- 
metric states equivalent under local unitary (LU) trans- 
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The main results of this paper are: a classification of 
the six local unitary stabilizer subalgebras (Lie algebras 
of the local unitary stabilizer subgroups) in Theorem [5] 
and for each of those algebras (except for the zero alge- 
bra), a classification of local unitarily distinct classes of 
states. In addition. Theorem [1] gives a structure theorem 
for stabilizer subalgebras of mixed states that generalizes 
a similar result for pure states in our earlier work [J]. 



The paper is organized as follows. We establish nota- 
tion and convention in Seetion|IIl Seetion lHll m'esents the 
identification of symmetric mixed states with polynomi- 
als in 3 variables and the basics of SO{3) representation 
theory (already given in our previous work [8[ but repro- 
duced here for the sake of self containment) needed for 
the analysis of stabilizers and their corresponding states. 
We give structure theorems for subalgebras of the local 
unitary algebra in general (Theorem [Ij , and stabilizer 
algebras of symmetric mixed states in particular (Theo- 
rem[2]), in Section ITVl Proofs of these theorems are given 
in the appendix. Five subsequent sections analyze and 
classify entanglement types of symmetric mixed states 
(local unitary equivalence classes) for five of the six stabi- 
lizer types (all but the zero algebra) given in Theorem [51 



II. PRELIMINARIES 

Let G ~ (S'[/(2))" denote the n-qubit local unitary 
group, where g ~ [91,92: ■ ■ ■ T9n) S G acts on a density 
matrix p by 

gpg'^ ■.= (9i®92®---® gn)p{gi «) 52 ® • • • «) 9,1)^. 

Let LG = ©"=1 su{2) denote the local unitary Lie alge- 
bra of infinitesimal transformations, where su{2) is the 
set of 2 X 2 skew-Hermitian matrices with trace zero. In 
order to study the actions of G and LG on density ma- 
trices, it is convenient to consider a larger vector space 
that contains the set of density matrices as a proper sub- 
set. Let W denote the 4-dimensionaI real vector space of 
2x2 Hermitian matrices. A convenient basis for W is 
{(To, cTi , (72 , era } , where ctq is the 2x2 identity matrix, and 
ci = cTx, (72 = (Jy, and (73 = (Tz are the Pauli matrices. 

The set of n-qubit density matrices is a proper subset 
of the vector space VV®", where every element p (whether 
or not p is positive or has trace 1) can be uniquely writ- 
ten in the form p = ^j s/C/, where / = 41*2 ■ ■ -in is a 
multiindex with ik = 0,1,2,3 for 1 < k < n, and aj 
denotes 

(jj = (jj^ (g) (7^2 <8) • • • (X" <yi„ , 

and the coefficients sj are real. When it is necessary to 
express a density matrix in terms of the computational 
basis, we use the notation p = J^i j Pi,J \I) {J\, where 
I ~ iii2 ■ ■ -in and J = J1J2 ■ ■ ■ jn are bit strings of length 
n. For / — iii2 ■ ■ ■ im we write Ik to denote the string 
*i*2 • • • *fc • • • *n obtained by complementing the kth in- 
dex, where i^ = ik + I (mod 2). Note that the same 
multiindex notation / = iii2 ■ ■ -in means ik = 0,1,2,3 
for the Pauli tensor expansion, and means ik ~ 0, 1 for 
the computational basis expansion, where context makes 
the meaning clear. 

An clement g <E SU{2) acts [2^ on at by 

(Ti h-^ ga^g^ (1) 

and the corresponding action of M in su{2) on at is 

at ^ [M, a,] = McTi - OiM. (2) 

An element 9 ~ (gi, 92, ■ ■ ■ , gn) in G acts on ct/ by 

ffCT/gt := (giai.gl) <E) ■ ■ ■ <^ (5«cri„gt ) (3) 

and M = (Mi, M2, . . . , M„) in LG acts on aj by 



[M,ai] 



fc=i 



(Til «) • • • «) [Mk^CTi^] 



(4) 



The local unitary stabilizer of p € W®" is the subgroup 
Stabp of G given by 

Stabp = {5 e G: gpg'' = p}. 



The corresponding local unitary stabilizer subalgebra, 
which we denote by Kp, is the Lie algebra of Stabp, given 
by 

Kp^{M eLG: [il/,p] = 0}. 

We will use the standard basis 

A = i(7z, B = i(7y, C = i(Tx 

for su{2). Given an element M £ su{2), we will write 
A/C") to denote the element (0, 0, . . . , Af, . . . , 0) in LG 
with M in the kth position and zero matrices in all other 
positions. 

We will make use of the following elementary 1-qubit 
Lie algebra actions on Pauli tensors. 



[A,a^.] = -2ay 

[A,(7y\ = 2(7^ 

[A,a,] = 



(5) 
(6) 
(7) 



Less trivial, but still elementary to check, are the follow- 
ing Lie algebra action calculations in standard coordi- 
nates 0]. 




= Y,ahJ)\i){j\ 
i,j 

= ^r;(/,J)|/)(J| 



(8) 
(9) 



I, J 



where C,{I, J), 77(1, J) are given by 

C(/,J) = 2ipi,jY,{-iy^c 



i 1.-^31 



liI,J) = i ^ickPi,,j-cepi,j,). 



k,e=i 

III. SYMMETRIC STATES AND 
POLYNOMIALS 

Given a permutation n of {l,2,...,n}, define P^^ to 
be the operator on Hilbert space that carries out the 
corresponding permutation of qubits. For example, 

P(23) IllOOO) = IIOIOO) 

^(23) (O-J «) 0-fe (g) (7;) P^23) == <^3 ® Ol ® (Tk 

Define a symmetrization operator on density matrices as 

Sym(p) = i-^P.pP-i. 
n! ^-^ 

■K 

An ?i-qubit symmetric density matrix p is one for which 
P.^pP~^ = p for all permutations tt, or equivalently, one 
for which Sym(p) = p. We denote by Sym" W the n-fold 



symmetric power of W defined in the previous seetion. 
It is the subspace of elements of W^" that are invari- 
ant under qubit permutation. Every n-qubit symmetric 
density matrix p is an element of Sym" W, and can be 
written 



^ XI ^"insns Sym {(J®'''' ® crl 



(10) 
where the sum is over non- negative integers ng, rii, n2, 
and 713 such that no + ni + n2 + na = n. The coefficients 
Crain2n3 ^.re real. The collection of n-qubit symmetric 
density matrices is a proper subset of Sym" W, since the 
latter contains Hermitian matrices that are not positive 
semi-definite, and Hermitian matrices for which the trace 
is not 1. 

Let M„[a;, 2/,z] be the set of polynomials of degree at 
most n in three variables a;, y, and z with real coeffi- 
cients. For each n, there is a linear map Fn : Sym" W -> 
M„[x, y, z] defined by 



^SymK- 



^0113 "l 



I— >■ X 



■"iy"^z"^ 



In this way, we may associate a polynomial of degree at 
most n with each n-qubit symmetric mixed state. The 
polynomial associated with (fTOj) is 



Fn{p) 






a;"i2/"=z"^ 



(11) 



For each n, the map Fn is an invertible linear map. 

Since Fn is a linear map, the polynomial for a mixture 
of symmetric mixed states is the mixture of the polyno- 
mials. 

FniPlPl +P2P2) ^PlFnipi) +P2Fn{P2) 

A product of polynomials R„[a:;,y, z] x IRm[2:,2/, -z] ^ 
M„+„i[x,y, 2] represents the symmetrized tensor product 
of states. 

Fn[pi)Fni{p2) = i^„+m(Sym(/9i ® P2)) 

Let g e SU{2). Define Tg : Sym" W ^ Sym" W to be 
the symmetric transformation of each qubit by g. 



Tg{p) ^ g^yig^f 



[x, y, z] 



i[x, y, z] to be the trans- 



(12) 



If we define Rg 

formation on polynomials defined by 

Rg{f)ix,y,z) ^ f {{x,y,z)^{g)) , 

where / is a polynomial, $ : SU{2) — ^ SO{3) is the ho- 
momorphism [2l| that associates a rotation in M'^ with 
each 2x2 unitary, and {x,y, z)^{g) denotes a row vec- 
tor multiplied by a 3 x 3 orthogonal matrix, then the 
following diagram commutes. 



Rn[x,y,z] 

-I 

Sym" W 



^n[x,y,z] 



Sym" W 



Homogeneous polynomials in three variables x, y, and 
z arc known to be reducible representations of 5*0(3) 
[131 H- If Vi is the irreducible representation of 50(3) 
with dimension 2/+1, then the homogeneous polynomials 
of degree p in three variables decompose into irreducible 
representations as 



Lp/2J 



2i- 



(13) 



IV. STABILIZER STRUCTURE FOR 
SYMMETRIC MIXED STATES 

We begin by considering an arbitrary subalgcbra 
K (not necessarily a stabilizer subalgebra) of LG = 
®r=i su(2). That is, K is a vector subspace of LG that is 
also closed under the Lie bracket. Theorem [T] below (the 
proof is given in the appendix) gives a structure theorem 
for how K decomposes as a direct sum of basic building 
blocks. It is convenient to give names to the following 
standard algebra types. 

Given a subset B C {1, 2, . . . , n} of qubit labels for a 
system of n qubits, let Ag denote the subalgebra 



Ab = 




M G su{2) 



of LG. In particular, note that Ag is isomorphic to su(2). 
Wc will call Ag the standard su{2) block subalgebra for 
the subsystem specified by B, and we will use the term 
su{2) block algebra for any algebra that is isomorphic to 
some Ae via local unitary group conjugation, that is, any 
algebra of the form 

(ffi «) .92 ® • • • «) 5„)^e(.9i «) ff2 «) • • • «) gn^ 

where g, G SU{2) for 1 < z < n. We will write iiTwcrnor 
to denote A{i,2,...,„} = {(Af,M, . . . ,M): M € sm(2)}, 
because it is the Lie algebra of the stabilizer group 
{{g,g, ■ ■ ■ ,g)- g & SU{2)} that defines the Werner states. 
We write A'p,oduct to denote the n-dimensional algebra 



K, 



product 



= I Y^ afc^C^) : ak G 



.fe=i 



and we write -ftToicko to denote the 1-dimensional algebra 

n 

that is the real span of the clement \J A^'^' . We use these 



k=l 



names because ^product is the stabilizer of the product 
state |00 • • • 0) (00 ■ ■ ■ 0|, and K^icke is the stabihzcr of the 
symmetric Dicke states Q . Finally, let /'S'ghz denote the 
(n — l)-dimensional algebra 



K 



GHZ 



n 

=1 



^a,ylW:a, eM,Xafe = 




which is isomorphic to m(I)*^"' ^\ The name of this al- 
gebra comes from the fact that it is the stabihzer of the 
GHZ stated. 

Theorem 1. Let K be a suhalgehra of LG = ®"^i su{2). 
The set {1, 2, . . . , n} decomposes as a disjoint union 

{1, 2, . . . , n} = 6i U 62 U • • • U 6p U 5 U 7^ 

such that the algebra K decomposes as a direct sum of 
algebras 



K 



with the following properties: 

(i) each Bi is an su{2) block subalgebra in qubits Bi, 

(ii) every element in the algebra S has zero coordinates 
in qubits outside of the qubit set S, 

(Hi) for every s in S , the set {Mg : ^j. Mf^ £ S} is a 
1-dimensional subalgebra of su{2), and 

(iv) every element in K has zero coordinates in qubits 
in TZ. 

Applying Theorem [T] to the stabilizer subalgebra Kp, 
where p is symmetric, leads to the following theorem, 
whose proof is in the appendix. 

Theorem 2. Let p be a symmetric mixed n-qubit state. 
There is an LU equivalent symmetric state p' such that 
the local unitary stabilizer subalgebra Kp> is one of the 
following. 

(a) all of LG 

(b) if Werner 
[CJ -tVpi'oduet 

(d) A'ghz 

(e) iiTDickc 

(f) the zero algebra 

In the following sections, we identify states that have 
the stabilizers given in Theorem [2 For stabilizer 
types (c), (d), and (e), we will make use of the following 
Lemma, whose proof is in the appendix. 

Lemma 1. Suppose that Kp is a stabilizer of type (c), 
(d), or (e), in Theorem\^ and that 

(gi®g2®---® gn)Kp{gi (g) 52 «> ■ ■ ■ ® ffn)''' = Kp 

for some gi,...,gn in SU{2). Then we have gkAgj, = 
±A for 1 < k < n. It follows that each g^ is either 
diagonal or antidiagonal. That is, gk is either of the 



form gk = 
real t. 







or gk 





At 



V. STABILIZER IS ALL OF LG 

The completely mixed state is stabilized by LG. In 
fact this is the unique state that has this stabilizer. One 
can see this as follows. 

Since C^ is an irreducible representation of 5*1/(2), it 
follows that (C^)*^" is an irreducible representation of 
G = (S'C/(2))" (see [H, Proposition 4.14). By Schur's 
lemma, any matrix that commutes with every g € G 
must be a scalar matrix ([1^, Theorem 1.10). 

Alternatively, it is straightforward to verify that if p 
has an off-diagonal element pu with / 7^ J, with say 
ik^jk, then (/|[A('=),p]|J) = {-iy''2ipij ^ (use ®). 
This rules out nonzero off-diagonal elements. Now sup- 
pose pii ^ Pikik foi' some k, where Ik is the same as 
/ in all positions but opposite in position k.. Then 
(/|[CW,p]|J) = i{pi,i, - pii) ^ (use ®). We con- 
elude that p is the completely mixed state. 



VI. STABILIZER IS Kwcmcr 

Let p be an n-qubit symmetric mixed state, and sup- 
pose that Kp = if Werner- It foUows that, on the group 
level, we have Stabp D {(g,g, ■■■,g)- g & SU{2)}. This 
is the class of Werner states. In [^ we prove that any 
symmetric Werner state must be of the form 

L«/2J 

p = ^ Cfe Sym (u®'^ ^ id®"-2fe j 



fe=0 



where 



1 CTr + cr„ (Ki cr„ + CTz 



for some real coefficients Ck.e- Further, we show that any 
two states with distinct choices of coefficients Ck,e are 
local unitarily inequivalent. 



VII. STABILIZER IS A'p^oduct 

It is clear that if p = ^^ sjcri where s/ = if ife = 
1,2 for every I < k < n, then Kp contains ifproduet- 
As long as there is some s/ 7^ for / 7^ (0,0,..., 0), 
we can apply (O to see that [C^''\ p] = 0, so we have 
Kp ~ A'product (the only possible algebra that contains 
ifproduet is all of LG). 

Conversely, suppose Kp = if product- Write p as 

P = 2Z ^-^^^ ® ^-^ '^ a2 ^■■"''1 ® ^'' 
J J 

+ X] ^^-^^2 ®(yj + ^ djas (g) a J 
J J 

where J is a multiindex of length n — 1. Then we have 



for some 







[^i'/'] =J2^-'2bja3 + 2cja2)®cFj. 



It follows that bj ~ cj — for all J. Thus we have 
established that Kp ~ A'pioduct if a-nd only if 

n 

p = ^CfeSym(af'=®Id®("-"') 

k=0 

for some real coefficients Cfc. 

In contrast to the case of symmetric Werner states 
in the previous section, LU equivalence classes of states 
whose stabilizer is -/^product rnay contain more than one 
state. For example, let p = Id/8 + aa^ ® <Jz ® (^z- Then 
p' = {iX)®^ p{~iX)®^ = Id/8 - aaz ®(Jz®(Jz is LU 
equivalent but not equal to p, and has the same stabi- 
lizer (here X denotes the Pauli matrix ax viewed as an 
element (7(2)). We claim that this type of sign change ac- 
counts for the entire LU class of a state for n > 3 qubits. 
That is, if p, p' are LU equivalent n > 3 qubit states with 
stabilizer Kp = Kpi ~ i^product, and p has correspond- 
ing polynomial Fn{p) ~ X]fc dkZ^ : then either p = p' or 
Fn{p') = J2ki~^)''^kz''- To see why this is true, suppose 
we have gpg'' = p' for some g = (gi, . . . ,g„) in SU{2)"'. 
We show in [3| that for n > 3, there exists an h € SU{2) 
such that p' = Th{p) = /i®"p(/it)®". By Lemma [U we 
have hAh'^ ~ ±A. Thus we have Rh{z) = ±z. 

Lemma [1] admits one additional possibility for 2-qubit 
states, and that is to have one of ^i, 172 be diagonal, and 
the other antidiagonal. An example is conjugation by 
ld®{iX) that takes p = ^-^-aa^i^a^ to p' — ^--aa^i^a^. 
It is easy to check that this exhausts all possibilities. 



Now suppose that p, p' are are LU equivalent n-qubit 
states of the form ([T^ for some n > 3. As was the case 
in the previous section, there is some g € SU{2) such 
that p' = Tg(p), and g is either diagonal or antidiago- 



nal. If g ~ 







, then conjugation by 5®" leaves 



p{n, k) fixed and performs a phase operation on po{a, /3) 
as follows. 



Tg{p{n,k)) 



p{n,k) 
PG(e'*" 



a, e 



'/3) 



If g = it r\ ' then Tg interchanges p{n, k) with 
p{n, n ~ k) and transforms pg{<^, P) as follows. 

Tg{p{n,k)) ^ p{n,n- k) 
Tg{po{a,f3)) = PG(e-^*"/3,e'*"a) 

The considerations of the preceding paragraph imply 
that if we exercise LU freedom to choose a, /3 real with 
a> (5, then (|14p is a unique representative of its LU class 
of states, with the only exception being the case when 
a ^ j3. In this case, p ~ J2k ^kpin, k) + dpQ{ly/2, l\/2) 
and p' ~ '^kc'uPin, k)+dpQ{l^/2, l-\/2) are interchanged 
by Tix- 

IX. STABILIZER IS /<Dicko 



VIII. STABILIZER IS Kgyiz 

First, we claim that if Kp = A'ghz, then pi.j = Q for 
J 7^ Z,/*^, where F denotes the bitwise complement of 
/. Indeed, \i J ^ I, F , then there exist qubit indices fc, I 
such that ife = jk and ii = j^. Then we have (apply (jS])) 

^ (/|[^W - A^'\p]\J) - C(/, J) - ±2*p/,j. 

Second, we claim that if /0/_/c ^ 0, then {I,F} = 
{00 • • • 0, 11 • • • 1}. Simply note that for any pair of posi- 
tions fe,^, we must have 

= (-1)'" -(-1)'* 

so ik ^ it- 

Let |V'(n,fc)) = II)®'' ® 10)®^"-'=) and let p{n,k) = 
Sym|V'(n,fc))(V(n,fc)|. Let |G(a,^)) = a|0)®"-h 
/3|1)®", and let PG(a,/3) = |G(a,/3)) (G(a,/3)|. The 
above two claims imply that if Kp = Kqyvl then p can 
be written as a mixture of the form 



P = ^ Ckp{n, k) + dpG{a, P) 

k=0 



(14) 



where Cq, ■ ■ ■ ,Cn,d are nonnegative real numbers that 
sum to 1. 



Let p be an n-qubit symmetric mixed state. We claim 
that g®"p(gt)®" = p for all diagonal g £ SU{2) if and 
only if Fn{p) is a linear combination of products of z'^ 
and {x'^ + y'^y, for < r + 2s < n. 

As for the case of the single su{2) block, we make use 
of the identification of symmetric states with polynomials 
given in section Hill 

If g is a diagonal element of SU{2), then ^{g) is 
a rotation about the z axis. Polynomials Fn{p) = 
'I2r s^r,sz'^{x'^ + if'Y (whcrc the sum is over nonnega- 
tive integers r, s such that r-\-2s<n) are invariant since 
z and x"^ + y^ arc invariant under rotations about the 
z axis. 

Conversely, suppose that g®^^p{g'^)®^ = p for all di- 
agonal g g SU(2). Then Fn{p) is invariant under ro- 
tations about the z axis. We now wish to decompose 
the homogeneous polynomials in three variables into a 
sum of irreducible representations of C/(l), and to note 
the dimension of the trivial representation, as it gives 
the space of polynomials invariant under rotations about 
the z axis. Each Vi in p^ decomposes into 21 + 1 one- 
dimensional irreducible representations of C/(l), of which 
one is the trivial representation. When expressed as a 
sum of irreducible representations of U{1), decomposition 
(fT3)) shows that the homogeneous polynomials of degree 
p in three variables contain [p/2j -I- 1 dimensions of the 
trivial representation. This is precisely the dimension of 
the space of homogeneous polynomials of degree p that 



can be produced by products of z^ and {x^ +2/^)" (*-e. 
the number of nonnegative integer pairs (r, s) for which 
r + 2s = p). The vector space M„[x, y, z] of polynomials 
of degree at most n is a direct sum of vector spaces of ho- 
mogeneous polynomials of degree p, for < p < n. Since 
we have accounted for the dimensions of each space of 
homogeneous polynomials, F„(p) must be a linear com- 
bination of the polynomials given. 

Now suppose that p, p' are LU equivalent, with 



Fnip) 



J2r,s^r,sZ''i^ 



V 



2\s 



and F„(p') 



J2r sK s^^ i^'^ 4- y'^y ■ For n > 3 qubits, the same 
analysis as for symmetric Werner states (Section IVI[) 
yields a diagonal or antidiagonal g £ SU{2) such that 



Tg{p) = p'. Since Rg{x^ + y^) 



y^ and Rg{z) = 



±z, the only nontrivial possibility is that Fn{p') = 
J2r si~'^y^r,sZ^{x'^ + y^Y ■ We Can discard the case 
n — 2, because the only possibilities for p are p = 
Id/4 -I- a{ax ®(Jx + (Jy®(Ty) and p — Id/4 -I- aa^ ® cTz ■ The 
first of these states has an su(2) block stabilizer, and the 
second has stabilizer -ftTproduct- 



X. SUMMARY AND CONCLUSION 

Table U shows a summary of LU classes of local uni- 
tary stabilizer algebra types and their corresponding LU 
classes of states. Having achieved LU classification for 
symmetric mixed states, it is natural to attempt further 
classes of mixed states. A natural avenue for investiga- 
tion is to take the stabilizer structures from Theorem [2] 
and try to classify their corresponding states, with the 
assumption of permutation invariance removed. For ex- 
ample, the case Werner states (stabilizer type (b)) would 
be of interest. 
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Appendix A: Proofs for Section IIVI 
1. Proof of Theorem [l] 

Let X be a subalgebra of LG — ®^i^isu{2). For 
each qubit i, let tt,; : LG — > su{2) denote the projection 
(Ml, Af2, . . . , Mi, . . . , Mn) — > Mi onto the ith direct sum- 
mand of LG. Given a set S of qubits, let ns denote the 
projection izg = ffises'i's onto summands in qubits S. 
Given M = (Mi, . . . , M„) in LG, let the weight of M, 
denoted wt(Af), be the number of i such that Mi is not 
zero. 

For each i £ {1,2, ...,n}, let m{i) be the minimum 
weight of all elements of K with nonzero weight in posi- 
tion i. For i such that m(i) > 0, let M(i) C K denote 



the set of elements of K that have nonzero projection in 
position i and have total weight m{i). That is, 

m{i) = min{wt(M): M e K : M, 7^ 0}, and 
M{i) = {Ad eK: M, ^ and wt(A/) = m{i)}. 



Proposition 1. Suppose diiRTTi{K) = 3, and let P,Q be 

elements of Mi 
only if Qj = . 



elements of M{i). For 1 < j < n we have Pj = if an 



Proof. Let U' be an element of K such that [//, Pi are 
independent elements of sw(2). Let V = P/\Pi\, let W = 
[U',V]/\[U',VU, and \etU=[V,W]. 

Observation 1. We have that Ui,Vi,Wi form an or- 
thonormal [2j1 basis of su{2) that satisfy [C/i,Fi] — Wi, 
[Vi,Wi] = Ui, and [Wi,C/i] = Vi. This is easy to see 
using the fact that su{2) is isomorphic to M.^, with the 
Lie bracket operation corresponding to the cross product, 
and the (rescaled) Hilbert-Schmidt norm corresponding 
to the standard norm on M.^. 

Observation 2. We have that U, V, W are elements of 
M{i). This follows from the observation that in K, we 
have wt([M,iV]) < wt(M) for ah M,N, so U,V,W ah 
have weight less than or equal to wt(P), and hence must 
have weight equal to wt(P) = m{i). Further, we see that 
[/, y, W must have zero and nonzero coordinates in the 
same positions as P. 

By Observation 1, we have that Qi is not a scalar 
multiple of at least one of Ui,Vi,Wi. Without loss 
of generality, suppose Qi,Ui are linearly independent. 
Let S = [Q,C/], so by construction, we have Si 7^ 0, 
S € M(i), and S has zero and nonzero coordinates in 
the same positions as U, and hence in the same posi- 
tions as P. If Qj 7^ 0, then Sj must be nonzero since 
wt(S') = wt{Q) = m{i) and S has the same zero and 
nonzero coordinate positions as P, which also has weight 
m,{i). Thus Uj 7^ 0, and therefore Pj 7^ 0. Conversely, 
if Pj 7^ 0, then Sj ^ 0, and therefore Qj 7^ 0. This 
concludes the proof. D 

Proposition 2. Suppose diimTi{K) — 3. Let B C 
{l,2,...,n} be the set of positions where all the ele- 
ments of M(i) have their nonzero coordinates, by virtue 
of Proposition [I) Then the following hold. 

(i) IfPeK and P, = for all j ^ B, then P € M{i) 
orP = 0. 

(a) Any such P lies in the linear span of U, V, W con- 
structed in the proof of Proposition [7J 

(m) The elements U, V, W satisfy [t/, V] = W, [V, W] = 
U, [W, U] =V so that the linear span of U, V, W is 
isomorphic to su{2). 

(iv) K = ttb{K) ® TTB-iK). 

(v) TTjs^K) is an su{2) block algebra. 



stabilizer 


LU Representative 




LU Nonuniqueness 


Pure State Example 


LG 


Id/2" 


unique 


none 


-f^ Werner 


Ln/2J 

^ Ck Sym Ko-^ ® Cr^r, + CTy (g) CTy + o-^ (g) (T^)®* (g) Id®<"" 
fc=0 

Cfc real, at least one Cfe 7^ with fc > 


■") 


unique 


singlet 


-^product 


^CfcSym(ar®Id®("-'=)) 

fc=0 

Cfc real, at least one Ck nonzero with A; > 


4 = (-i)'cfe 


product state 


^GHZ 


n 

'^Ckpin,k) + dpGG{a, P) 

fe=0 

a> P>0 


c'fe = C„_fc 


GHZ 


-K'Dickc 


J2 ^= Sym ((a, ® a, + a^ ® a,)®= ® af ® Id®"-"- 

2r + s<n 

6r,s real, at least one b^.s / with r > 


-") 


&;.s = i-iy-br.s 


Dicke state 



TABLE I: Summary of LU classes of stabilizers and states. Each of the three last stabilizer types has precisely two LU 
inequivalent states of the form given in the column 'LU Representative'. The column 'LU Nonuniqueness' gives conditions on 
the coefficients for the alternative state. 



cW. 
—aUi 



J, and let Q 



Proof, (i). If P = or if P; 7^ we are done, so suppose 
Pi = 0, and let j 7^ i be an clement of B such that 
Pj 7^ 0. Let C/, V^ W be constructed as in the proof of 
Proposition [T] Since Uj,Vj,Wj are nonzero, they are 
independent elements of su{2). Write Pj = aUj + hVj + 
P - aU ~hV - cW. Then Q, = 
bVi - cWi 7^ 0, so Q e M{i). But Qj = 0, and 
this contradicts Proposition [T] We conclude that Pi = 
is impossible, so P S M{i). 

(ii). Using (i), for P 7^ we can write Pi — aUi + bVi + 
cW, and \ct Q = P - aU - bV - cW. Then Qi = 0, so 
Q = 0, again by (i). Thus P ^ aU + bV + cW. 

(iii). We only need to show that [Uk,Vk] = Wk, 
[Vk,Wk] = Uk, and [Wk,Uk] = Vk for k e B. By con- 
struction, we have [V, W] = U, and therefore we have 
[Vk,Wk] = Uk. We also have [Wk,Uk] = UkVk and 
[C^A;,V/c] ~ PkWk for some scalars ak,/3k- Let N ~ 
[V, [U, V]], and let M = N-U. Then Nk = [Vk,/3kWk] - 
PkUk- Since Mi = 0, we must have Mk = [Pk - l)C//c == 
by (i), so we must have /3fc = 1 for all /c G S. A similar 
argument shows that a^ = 1 for all k G B. 

(iv). Let M G A, and write M = P + N, where P = 
TTBiM) and N = M - P. We will show that P G ^. We 
have [M, X] = [P, X] for X = U,V, W, so we can use (ii) 
to write [P, X] as a linear combination of U, V, W. 

[P,U] = auU + buV + cuW 
[P,U] = avU + bvV + cvW 
[P,W] = awU + bwV + cwW 



coefficients in the second set of equations yields 



For k G B,we can write Pk 
the relations (iii), we have 

[PkMk] = 
[Pk,Vk] 
[Pk,Wk] 



^akUk + bkVk + CkWk- Using 

-bkWk + CkVk 
dkWk - CkUk 
-akVk +bkUk- 



ak = 


= -bw = cv 


bk = 


- ^cii = aw 


Ck = 


-- = —ay = 6(7 



Since this holds for all k, we have P = cyU+awV+buW, 
and so P G X, as desired. 

(v). Let b e B. Because Ub,Vb,Wb satisfy [£/&, Vf,] = 
Wt, [Vb,Wb] = Ub, and [Wb,Ub] = Ub, it must be that 
Ub,Vb,Wb form an orthonormal basis for su{2) (again 
use the fact that su{2) with its Lie bracket and Hilbert- 
Schmidt norm is isometrically identified with M.^ with the 
cross product and standard euclidean norm). Because 
the adjoint representation SU{2) — >• SO{su{2)) is surjec- 
tive, we may choose hb G SU{2) such that hbUbh\ = A. 
Now choose a real number tb such that e**''V{,e~'*'' = B. 
Now let gb = e**''/!^. For b' ^ B, let gb' = Id. We now 
have 



{®^g^)U{®^gi)^ 

{.®igi)V{(^igi)'^ 
{®igi)W{®^gi)'^ 



{A,A,...,A) 
{B,B,...,B) 
(C,C,...,C) 



Equating fcth coefficients of the first set of equations with 



and TTig{K) has been "aligned" with Ag, as desired. D 

To complete the proof of Theorem [l] for each qubit 
label i for which dim7ri(iir) = 3, let B C {1, 2, . . . , n} be 
the set of positions where all the elements of M{i) have 
their nonzero coordinates, as in Propositions [T] and [51 
Then tt]3(K) is an su(2) block summand for K by part (v) 
of Proposition!!] For any qubit j outside of the qubit sets 
for su{2) blocks, we must have dim7rj(A') — 0,1, for if 
there are two independent vectors in ttj(K), then there 
must also be a third coming from the bracket of the two 
independent elements, and so j would be a qubit for an 
su{2) block. We define S to be the set of qubits j for 
which dim7rj(i^) = 1 and define TZ to be the remaining 



qubits £, where we have dim7r£(ii') = 0. With these 
definitions, we clearly have the desired decomposition of 
Theorem [TJ 



2. Proof of Theorem [2] 



but not through an LU transformation applied uniformly 
to each qubit. In fact, the singlet is invariant under any 
LU transformation applied uniformly to each qubit. 

Suppose p is a 2-qubit symmetric mixed state whose 
stabilizer subalgebra is an su{2) block subalgebra. Writ- 
ing 




Let p be a symmetric mixed state, and let Kp be its 
local unitary stabilizer subalgebra. By Theorem [H we 
can decompose Kp as a direct sum 



(Al) 



of su{2) blocks Bi and an algebra S where projections 
into each qubit summand are 1-dimensional. 

To begin, note that permutation invariance implies 
that the dimension of the projection of Kp into any su{2) 
summand of LG must be the same for all qubits. Thus 
we can have only the zero algebra (possibility (f) of The- 
orem [2]), an algebra S, or a sum of su{2) blocks. 

Next, we claim that there are only two possibilities for 
a sum of su{2) blocks. One extreme is to have n su(2) 
blocks, each in 1 qubit (possibility (a) of Theorem [2|). 
The other is to have a single su{2) block in all n qubits. 
We rule out the intermediate possibilities, that is, having 
two or more su{2) blocks, at least one of which involves 
two or more qubits, as follows. Suppose Bi contains 
qubits i and j, with corresponding su{2) block algebra 
J3i, and B2 contains qubit k ^ i,j, with corresponding 
su(2) block algebra B2. Transposing qubits j, k does not 
affect p, so Kp also has an su{2) block algebra B3 that 
contains qubits i, k. But then the qubit sets for Bi and 
B^ both contain i. This contradicts the fact that qubit 
sets for su{2) blocks arc disjoint, as shown in Proposi- 
tion H 

To complete part (b) of Theorem [5J we consider two 
cases. 

Part (b) case (i) 71 > 3. We claim that if Kp is a single 
su{2) block in all n > 3 qubits, then in fact Kp is the 
standard su{2) block algebra Am 2,....n}- Suppose on the 
contrary that there is an element M ~ ^^ M^ in Kp 

with A'U 7^ Mj. Let M' = Nlf + M^^^ + T.k^^. ^f'^ be 
the element obtained from M by transposing the i, j co- 
ordinates. Then M-M' = {M, - Mj)'-'^ + {Mj - M,)^^'^ 
is also in Kp, but this element has weight 2 < n. This 
contradicts Proposition [2] that says all elements in an 
su{2) block must have full weight n. We conclude that a 
single su{2) block stabilizer for a symmetric mixed state 
of n > 3 qubits must be the standard su{2) block algebra. 
Part (h) case (ii) n = 2. We claim that any two-qubit 
symmetric mixed state with an su(2)-block stabilizer is 
LU equivalent to another symmetric state with a stan- 
dard SM(2)-block stabilizer. This LU equivalence may 
not be achievable through the same unitary operation on 
each qubit. For example, the pure state |01) -I- |10) is a 
symmetric state with a non-standard sw(2)-block stabi- 
lizer. It is LU equivalent to the singlet state |01) — |10), 



P — / J Sij<Ti Yy uj, 
ij=0 



) ai 



we see that 



sio = S20 = S30 = soi = S02 = so3 = 

because the three components (sio, S20, S30) transform 
like a 3-dimensional real vector under the rotations of 
3-dimensional space produced by the first qubit in the 
su{2) block subalgebra. Since p is invariant under the 
action of this su{2) block subalgebra, it must be that 
(siOi S207 Sao) = 0. A similar argument holds for the sec- 
ond qubit. By performing singular value decomposition 
on the real 3x3 matrix of elements Si-^i^ with ii and 12 
ranging over {1,2,3}, we can find a local unitary trans- 
formation such that p' = (pi (X) g2)p{gi ^ 32)^ has the 
form 

p' = - (o-Q ® ctq + acTi ®ai+ ba2 ® a2 + ca^ ® 0-3) . 

Let 

M = aiA(i)-f/3iB(^)+7iC*^^+a2A(^)-h/32B*^^+72C(2). 

This M stabilizes p' if and only if 

aic ~ a2b 
aib = a2C 

Pia = P2C 
/3ic = /32a 
716 = 72a 
71 a = 72^. 

Since there is, by assumption, an element M with ai, /3i, 
71, a2, h, and 72 all nonzero, we must have 

\a\ = \b\ = \c\. 

If a, 6, and c are not identical, then two have the same 
sign and the last has the opposite sign, and we can do 
a unitary transformation on one qubit that produces a 
■n rotation about the x, y, or z axis, resulting in an 
LU equivalent state that is symmetric with the standard 
su{2) block subalgebra as its stabilizer. This completes 
part (b), case (ii). 

Now we suppose that Kp = S, so that the projection 
of Kp into each LG summand su(2) is 1-dimensional. 
We wish to show that one of possibilities (c)-(e) of the 
statement of the Theorem holds. We do this by cases. 



fore we have M - 
in Kp. Let N = 



If Kp contains a weight one element M^^\ then by 
permutation invariance, Kp must contain M^^' for 1 < 
k < n. This estabUshes possibiUty (c). 

Now suppose that Kp contains no weight one elements. 
There are two cases: there may exist a nonzero element 

M = (Mi,M2,...,A//„) = Efe^4''^ ii^ Kp such that 
Mi y^ Mj for some i,j, or there may not. Suppose the 
first case holds. For any pair i,j of qubit positions, we 
have the element M' = il/j'^ + il/f^' + Efc^^j ^4''^ ob- 
tained by transposing coordinates i,j, also in Kp. There- 

M' = (M, - Mj)^*^ + {Mj - Mif''' also 
Mi — Mj . We have n — 1 linearly in- 
dependent weight two elements N^^' — N^'^' , 2 < k < n, 
so the dimension of Kp is at least n — 1. If the dimen- 
sion of Kp is greater than n — 1, then there is an element 
in Kp of the form M" = X]fc=i o-kN^''^ that docs not 
lie in the span of the elements iY*^') — A^'^-'^. But then 
we would have the weight one element (X]fe=i '^fe) ^'"'^^ ~ 
M" + J2l=2°'i {^'-^^ - ^'''') in Kp, which contradicts 
our assumption. We conclude that Kp is the real span 
of {7V(i) - A^(*)}"^2- We may take N to have norm 1, 
and we may choose g G SU{2) such that gNg^ = A. 
Let p' = ,g®"p(5®")t. Then Kp, is the real span of 
{A^') - ^(•'^}i<ij<„, which is the stabilizer (d) of Theo- 



Mj for all 

= S. Let 

be a nonzero element in Kp. Nor- 



rem[2l 

The last case to consider is where Mi - 
i,j, for all M = (Mi, M2, . . . , Af„) in Kp 
M = {N,...,N) 
malize N and choose g E SU{2) to diagonalize N so 
we have gNg"^ = A. Then the state p' = ^'^"^(g'f)®" 
has 1-dimensional stabilizer that is the real span of 
^(1) + A(2) + h A("). This is type (e) in Theorem H 

This concludes the proof of Theorem [2l 



3. Proof of Lemma [T1 

Suppose that Kp is a stabilizer of type (c), (d), or (e), 
in Theorem [21 and that 



(51 (g) 52 ® ■ ■ ■ <8i 9n)Kp{gi ® 52 fgi ■ ■ • ® ffn) 



for some gi, . . . ,gn in SU{2). For M in any of the three 
stabilizer types, the assumption that {®kgk)M{®k{gk) G 
Kp implies that gkAg\ ~ ±A for each k. It is easy to 

a —h 

b a 

then either a or b must be zero. Thus we conclude that 
g is either diagonal or antidiagonal, as desired. 



check directly that if g = 



satisfies gAg^ = ±A, 
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When we identify the element aa^ + bay + ca^ in W with 
i{aax + bay + ca^) in su{2), equation ((T} is the adjoint 
action of 5(7(2) on its Lie algebra su{2), plus a trivial 
action on the real linear span of ctq. Equation ^ is the 
corresponding adjoint action of the Lie algebra su(2) on 
itself, plus again the trivial action on ao. Equations Q 
and @ are the natural extensions of ([TJ and ^ to tensor 
products. See, for example, [19l |. 

"I> is given in a natural way by the adjoint action SU{2) — >■ 
SO{su{2)), so that $(p)(M) = gMg\ and we iden- 
tify su{2) with R^ by yl ^ (1,0,0), B ^ (0,1,0), 
C -^ (0,0,1). See ^. 

To be precise, the cited work considers the representa- 
tion C[x,y,z]. In the case of 5*0(3), the irreducible sub- 
modules of this complex representation are in one-to-one 
correspondence with the real irreducible submodules of 
R[a:,2/,2:] via complexification. See [l9l . Ch.2 Sec. 6]. 
We use a rescaled Hilbert-Schmidt norm [Mj = 



10 

(2tr(MtM))^/2 for an element M of su{2). 



